Bose-Einstein Condensation of Dark Matter Axions 
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We show that cold dark matter axions thermalize and form a Bose-Einstein condensate. We 
obtain the axion state in a homogeneous and isotropic universe, and derive the equations governing 
small axion perturbations. Because they form a BEC, axions differ from ordinary cold dark matter 
in the non-linear regime of structure formation and upon entering the horizon. Axion BEC provides 
a mechanism for the production of net overall rotation in dark matter halos, and for the alignment 
of cosmic microwave anisotropy multipoles. 
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Several authors have proposed that the dark matter of 
the universe is a Bose-Einstein condensate (BEC) [l|, Q • 
The axion is sometimes mentioned in this context. In- 
deed the axion is a boson and a cold dark matter can- 
didate, and cold dark matter axions are known to have 
huge phase space density. But, as far as we are aware, 
it has never been shown that dark matter axions form a 
BEC. Their phase space density is certainly large enough 
but they will only form a BEC if they reach thermal equi- 
librium. This may see unlikely because the axion is very 
weakly coupled. Below we find that dark matter axions 
do form a BEC, marginally because of their self inter- 
actions but certainly as a result of their gravitational 
interactions. No special assumptions are required. 

Shortly after the Standard Model of elementarypar- 
ticles was established, the axion was postulated [3j] to 
explain why the strong interactions conserve the discrete 
symmetries P and CP. For the purposes of this paper the 
action density for the axion field <p(x) may be taken to 
be 
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where m is the axion mass. The self-coupling strength is 
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in terms of the axion decay constant / and the masses 
m u and m,d of the up and down quarks. In Eq. {U, tne 
dots represent higher order axion self-interactions and 
interactions of the axion with other particles. All axion 
couplings and the axion mass 
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are inversely proportional to /. / was first thought to be 
of order the electroweak scale, but its value is in fact arbi- 
trary [4]. However, the combined limits from unsuccess- 
ful searches in particle and nuclear physics experiments 
and from stellar evolution require / > 3 • 10 9 GeV 

Furthermore, an upper limit / < 10 12 GeV is provided 
by cosmology because light axions are abundantly pro- 
duced during the QCD phase transition In spite of 



their very small mass, these axions are a form of cold 
dark matter. Indeed, their average momentum at the 
QCD epoch is not of order the temperature (GeV) but 
of order the Hubble expansion rate (3 • 10~ 9 eV) then. In 
case inflation occurs after the Peccei-Quinn phase tran- 
sition their average momentum is even smaller because 
the axion field gets homogenized during inflation. For a 
detailed discussion see ref. [7[. In addition to this cold 
axion population, there is a thermal axion population 
with average momentum of order the temperature. 

The non-perturbative QCD effects that give the ax- 
ion its mass turn on at a temperature of order 1 GeV. 
The critical time, defined by m(ti)ti — 1, is ti ~ 
2 ■ 1CP 7 sec (//10 12 GeV) 3. Cold axions are the quanta 
of oscillation of the axion field that result from the turn 
on of the axion mass. They have number density 
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where a(t) is the cosmological scale factor. Because the 
axion momenta are of order j- at time t\ and vary with 
time as a(t) -1 , the velocity dispersion of cold axions is 



mti a(t) 



(5) 



if each axion remains in whatever state it is in, i.e. if 
axion interactions are negligible. Let us refer to this case 
as the limit of decoupled cold axions. If decoupled, the 
average state occupation number of cold axions is 
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Clearly, the effective temperature of cold axions is much 
smaller than the critical temperature 
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for BEC. Axion number violating processes, such as their 
decay to two photons, occur only on time scales vastly 
longer than the age of the universe. The only condition 
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for axion BEC that is not manifestly satisfied is thermal 
equilibrium. 

Axions are in thermal equilibrium if their relaxation 
rate T is large compared to the Hubble expansion rate 
H(t) = ir. At low phase space densities, the relaxation 
rate is of order the particle interaction rate T s — naSv 
where a is the scattering cross-section. The cross-section 
for ip + *p — > <p + ip scattering due to axion self interaction 
is in vacuum 

a = ^^~1.5-lQ- W5 cm 2 (-^-^Y . (8) 
64tt m 2 VIO- 5 eVJ w 

If one substitutes <7o for cr, T s is found much smaller than 
the Hubble rate, by many orders of magnitude. However, 
in the cold axion fluid background, the scattering rate is 
enhanced by the average quantum state occupation num- 
ber of both final state axions, a ~ (TqN 2 '', because energy 
conservation forces the final state axions to be in highly 
occupied states if the initial axions are in highly occupied 
states. In that case, the relaxation rate is multiplied by 
one factor of M [|| 

r - n a Sv JV . (9) 

Combining Eqs. (I4ll6l8l). one finds r(ti)/H(ti) ~ 0(1), 
suggesting that cold axions thcrmalize at time t\ through 
their self interactions, but only barely so. 

It may seem surprising that the huge and tiny fac- 
tors on the RHS of Eq. §§§ cancel each other. In fact 
the cancellation is not an accident. Consider a generic 
axion-like particle (ALP) whose mass m and decay con- 
stant / are unrelated to each other. Its self interaction 
2 

coupling strength A ~ -jr. Cold ALPs appear at a time 
t\ ~ ~ with number density n{t\) ~ f 2 m, and velocity 
dispersion Sv(t±) ~ 1. Substituting these estimates in 
Eqs. ([6]), (JHJ) and ((9]), one finds that the thermalization 
rate is of order the Hubble rate at t\, for all / and m. 

A critical aspect of axion BEC phenomenology is 
whether the BEC continues to thermalize after it has 
formed. Axion BEC means that (almost) all axions go to 
one state. However, only if the BEC continually rether- 
malizes does the axion state track the lowest energy state. 

The particle kinetic equations that yield Eq. ((9]) are 
valid only when the energy dispersion ^m(Sv) 2 is larger 
than the thermalization rate [8J. After t\ this condition 
is no longer satisfied. One enters then a regime where 
the relaxation rate due to self interactions is of order 
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T\(t)/H(t) is of order one at time t\ but decreases as 
t a(t)~ 3 afterwards. Hence, self interactions are insuffi- 
cient to cause axion BEC to rethermalize after t\ even if 
they cause axion BEC at t\. However gravitational in- 
teractions, which are long range, come in to play. The 
relaxation rate due to gravitational interactions is of or- 
der 



where £ ~ (mSv) 1 is the correlation length. T g (t)/H(t) 
is of order 4 • 10™ 8 (//10 12 GeV)§ at time t\ but grows as 
ia _1 (t) oc a(t). Thus gravitational interactions cause the 
axions to thermalize and form a BEC when the photon 
temperature is of order 100 eV (//10 12 GeV)^. 

The process of axion Bosc-Einstein condensation is 
constrained by causality. We expect overlapping conden- 
sate patches with typical size of order the horizon. As 
time goes on, say from t to 2t, the axions in i-size conden- 
sate patches rethermalize into 2i-size patches. The corre- 
lation length is then of order the horizon at all times, im- 
plying Sv ~ instead of Eq. (J5J, and T g /H cx t 3 a~ 3 (t) 
after the BEC has formed. Therefore gravitational inter- 
actions rethermalize the axion BEC on ever shorter time 
scales compared to the age of the universe. 

We now consider what implications axion BEC has for 
observation. The axion field may be expanded in modes 
labeled a: 
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where the 3>s(x) are the positive frequency c- number so- 
lutions of the Heisenberg equation of motion for the axion 
field 

D»Dp(p(x)=g^[dpd v -r^dxMx)=m 2 tp(x) . (I 3 ) 

and the as and at are creation and annihilation operators 
satisfying canonical commutation relations. We neglect 
the self- interaction term — gAc/? 3 , which would otherwise 
appear on the RHS of Eq. (flU|) , because it is of order 
-pip, where p is the axion density, and hence smaller by 

the factor ^^pj-^ than the relevant terms (of order 

Y<p) m that equation. When the self- interactions are in- 
cluded, one finds an instability in the axion BEC towards 
the formation of droplets. The analog of the sound speed 
is imaginary in this case because the self interaction 
is attractive. However, the rate of droplet formation is 
negligibly small compared to the Hubble rate. The gravi- 
tational forces always dominate over the self interactions 
except briefly after the cold axions first appear at time 

h 

Except for a tiny fraction, all cold axions go to a single 
state which we label a = 0. The corresponding &o(x) 
is the axion wavefunction. In the spatially flat, homoge- 
neous and isotropic Robertson- Walker space-time, 

$ = -A3- e~ mt (14) 
a(t) 2 

where A is a constant. The state of the axion field is 
\N >= (1/VM) (4)^10 > where |0 > is the empty 
state, defined by as |0 > = for all a, and N is the 
number of axions. The expectation value of the stress- 
energy-momentum tensor is 

<N\T MV \N> = N[d^* d„$ 
+d v $* d fl $ + gfW (-d x %d x ^ -m 2 ^ )} (15) 
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Again we neglect the self-interaction term. 

Consider first the behavior of axion BEC in a flat 
Minkowski space-time. Since the axions are non- 
relativistic, &o(x) — e~ lrntl $>(x) with 4"(a;) slowly vary- 
ing. Neglecting terms of order ^ dt compared to terms 
of order one, Eq. (TIB"]) becomes the Schrodinger equation: 
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It is useful [Io| to write the wavefunction as 
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In terms of B(x, t) and /3(x, t) the energy and momentum 
densities are (J, k = 1,2,3) Too = P = m(B(x,t)) and 
Toj = — pvj = — (B(x, t)) 2 dj/3, in the non-relativistic 
limit. The velocity field is therefore v{x,t) — — V(3(x,t) 
fTpj ] . Eq. (fTB"|) implies the continuity equation and the 
equation of motion 



d t v k +v j djv k = -V<j 
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Following the motion, the stress tensor is 
Tj k = pvjV k + -^^{-djpdkp - 5 3 k^ 2 p) 
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For ordinary cold dark matter (CDM) the last terms on 
the RHS of Eqs. HID and (20]) are absent. 

To compare axion BEC with CDM we divide the obser- 
vations into three arenas: 1) the behaviour of density per- 
turbations on the scale of the horizon, 2) their behaviour 
during the linear regime of evolution within the horizon, 
and 3) their behaviour during the non-linear regime. We 
first discuss arena 2 where CDM provides a very suc- 
cessful description. Neglecting second order terms, the 
perturbation in the stress tensor implied by Eq. ([20)) is 



ST jk = -Sjk^-^-V 2 d(x,t) 
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where po(t) is the unperturbed axion density and 
S(x, t) = Because the RHS of Eq. fH]) is propor- 

tional to the Kronecker symbol and the RHS of Eq. (TT8|) 
is a gradient, vector and tensor perturbations are not af- 
fected by the additional forces associated with the axion 
BEC. Only the scalar perturbations are affected. The 
scalar perturbations are conveniently described in con- 
formal Newtonian gauge [IH where the metric is 



Conservation of energy and momentum in this back- 
ground implies the first order equations 
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where H — -j^jf • The equations for CDM are recovered 
by letting m — > oo. The RHS of Einstein's equations are 
modified by the addition of STjk to the stress tensor, but 
this modification does not play a role in our discussion 
because it is suppressed, relative to the leading terms, by 

the factor (^^-j , where fc p h is the physical wavevector 
of the perturbation. 

It is clear from Eqs. (|2"3")) that axion BEC differs from 
CDM on small scales only. For scales that are well within 
the horizon (fc p h >> H), Eqs. (|23|) plus Einstein's equa- 
tions imply 



d 2 S + 2Hd t 5 - 
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for the Fourier components S(k,t) of 5(x,t). k = ak p h 
is co-moving wavevector. We assumed <p = —ip which 
is almost always the case [l~H and certainly valid during 
the matter dominated era. Eq. (|24|) shows that the axion 
BEC has Jeans length 

kj 1 = {16nGpm 2 )-i 

^ l-02-10 14 cn 1 (^) A ( 10 " 29 ; /Cm3 )(25) 

The Jeans length is small compared to the smallest scales 
(~ 100 kpc) for which we have observations on the behav- 
ior of density perturbations in the linear regime. Thus 
axion BEC and CDM are indistinguishable in arena 2 on 
all scales of observational interest. 

In the non-linear regime of structure formation (arena 
3) and in the absence of rethermalization, the relevant 
equations are 



d t p + V ■ [pv) = 
d t v + (v ■ V)w 



V x v = 
-V-0 - v<? 



(26) 



-(l+2ip(x,t))dt 2 +a(t) 2 (l+2cf)(x 1 t))dx-dx . (22) 



Eqs. (|26p are equivalent to the Schrodinger equation for 
particles in a Newtonian gravitational field. Axion BEC 
and CDM differ in that the — Vq term is absent from the 
force law for CDM. However, as was shown by numerical 
simulation . and as is expected from the WKB ap- 
proximation, the differences occur only on length scales 
smaller than the de Broglie wavelength. Since the ax- 
ion de Broglie wavelength (of order 10 meters in galactic 
halos) is negligbly small compared to all length scales of 
observational interest, we again find that axion BEC and 
CDM are indistinguishable when there is no rethermal- 
ization of the BEC. 
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However, we found that gravitational interactions do 
rethermalize the axion BEC continually so that the axion 
state tracks the lowest energy state. This is relevant to 
the angular momentum distribution of dark matter ax- 
ions in galactic halos. The angular momentum of galaxies 
is caused by the gravitational torque of nearby galaxies 
early on when protogalaxies are still close to one another 
fl3j ]. The angular momentum distribution acquired by 
the dark matter particles determines the structure of the 
inner caustics that they form in galactic halos 14 , ll| . If 
that distribution is characterized by net overall rotation, 
implying Vx«^0, the inner caustics are a set of "tricusp 
rings" [l4J . If the velocity field is irrotational (V x v = 0) , 
the inner caustics have a tent-like structure [l5( quite dis- 
tinct from that of tricusp rings. Evidence has been found 
for tricusp rings [l6| , as opposed to the tent-like caustics 
of the V x v = case. This raises a puzzle for CDM. In- 
deed one can show [3] that the velocity field of ordinary 
cold dark matter, such as WIMPs, remains irrotational 
as it is the result of gravitational forces which are propor- 
tional to the gradient of the Newtonian potential. The 
puzzle is solved if the dark matter is an axion BEC which 
rethcrmalizcs while tidal torque is applied to it. Indeed, 
the lowest energy state for given total angular momen- 
tum is one in which each particle carries an equal amount 
of angular momentum. In that case there is net overall 
rotation. V x v ^ is accomodated in the BEC through 
the appearance of vortices. The phenomenon is observed 
in quantum liquids and well understood [10( ■ 

Finally we consider the behaviour of density pertur- 
bations as they enter the horizon (arena 1). Here too 
axion BEC may differ from CDM. The CDM perturba- 
tions evolve linearly at all times. The axion BEC pertur- 
bations do not evolve linearly when they enter the hori- 
zon because the condensates which prevailed in neighbor- 
ing horizon volumes rearrange themselves, through their 
gravitational interactions, into a new condensate for the 
expanded horizon volume. This produces local correla- 
tions between modes of different wavevector since the 
perturbation of wavevector k, upon entering the horizon, 
is determined by the perturbations of wavevector say \k 
in its neighborhood. We propose this as a mechanism 
for the alignment of CMBR anisotropy multipoles 17 1 
through the integrated Sachs- Wolfe (ISW) effect. Un- 
like CDM, the ISW effect is large in axion BEC because 
the Newtonian potential %j) changes entirely after enter- 
ing the horizon in response to the rearrangement of the 
axion BEC. 

We conclude that a case can be made that a large frac- 
tion of the dark matter is axions. Although the QCD ax- 
ion is best motivated, a large class of axion-like particles 
has the properties described here. 
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